To appear in ApJL 



Dusty disks at the bottom of the IMF 

Alexander Scholz 

SUPA, School of Physics & Astronomy, University of St. Andrews, North Haugh, St. 

Andrews, KYI 6 9SS, United Kingdom 

asllOQst-andrews .ac.uk 

Ray Jayawardhana 

Department of Astronomy & Astrophysics, University of Toronto, 50 St. George Street, 

Toronto, ON M5S 3H4, Canada 

rayj ayOastro . utoronto . ca 
ABSTRACT 

'Isolated planetary mass objects' (IPMOs) have masses close to or below the 
Deuterium-burning mass limit (~ 15 Mjup) - at the bottom of the stellar ini- 
tial mass function. We present an exploratory survey for disks in this mass 
regime, based on a dedicated observing campaign with the Spitzer Space Tele- 
scope. Our targets include the full sample of spectroscopically confirmed IPMOs 
in the a Orionis cluster, a total of 18 sources. In the mass range 8 ... 20 Mjup, we 
identify 4 objects with > 3a colour excess at a wavelength of 8.0 /um, interpreted 
as emission from dusty disks. We thus establish that a substantial fraction of 
IPMOs harbour disks with lifetimes of at least 2-4 Myr (the likely age of the clus- 
ter), indicating an origin from core collapse and fragmentation processes. The 
disk frequency in the IPMO sample is 29±}3% at 8.0 /jm, very similar to what has 
been found for stars and brown dwarfs (~ 30%). The object SOri70, a candidate 
3 Mjup object in this cluster, shows IRAC colours in excess of the typical values 
for field T dwarfs (on a 2a level), possibly due to disk emission or low gravity. 
This is a new indication for youth and thus an extremely low mass for SOri 70. 

Subject headings: stars: circumstellar matter, formation, low-mass, brown dwarfs 
- planetary systems 
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Introduction 



Over the past five years, it has been firmly established that the evolution of brown 
dwarfs from 1-10 Myr follows the blueprint known from T Tauri stars: They show evi dence 
for ongoing gas accretion and outflows (e.g. iMuzeroUe et al.ll2005l : iMohantv et al.ll2005), and 



Jayawardhana et al 



2003 



harbour circum-su b-stellar disks with lifetimes of 5-10 Myr (e.g 
Scholz et al.l 120071 ). not vastly different from stars. This has often been interpreted as an 
indication for a common origin of stars and brown dwarfs - sub-stellar objects as the nat- 
ural extension of the stellar initial mass function (IM F) and the lowest-rn ass outcome of 
core collapse and fragmentation processes (see review by iLuhman et al.ll2007l . and references 
therein). At the same time, ongoing deep surveys have pushed the sensitivity limits down to 
masses below the Deuterium burning limit at ~ 15 Mjnp, finding evidence fo r a cluster popu- 
lation with masses comparable to giant planets (IZapatero Osorio et al.ll2000l : iLucas fc Roche 
2000l ). These objects at the bottom of the IMF are inconsistently called planemos, sub-brown 
dwarfs, cluster planets, or isolated planetary mass objects in the literature; in the absence 
of a definite nomenclature, we will use the latter term (abbreviated IPMO) in the following. 

The nature and origin of IPMOs have been debated intensely in recent years. Most 
of the problems raised in the discussion s about the formation of brown dwarfs apply even 
more acutely to IPMOs (see review by IWhitworth et al.l 120071 ). For example, it is under 
debate whether IPMOs can form directly from the collapse of ultra-low mass cores; instead, 
they may be giant planets ejected from circumstellar disks or stellar embryos ejected from 
mini-clusters or decaying multiple systems. 

As of today, the observational constraints on the nature of IPMOs are very limited. Ev- 
idence for accretion or disks has been reported for a few young sources with masses around 



the Deuterium burning l imit (Natta et al. 20021: [zapatero Osorio et al. 2002a; Luhman et al. 



20051 : lAUers et al.l l2006l : IJayawardhana fc Ivanovl l2006l ) . The next step is clearly to estab- 
lish frequency and lifetimes of IPMO disks, and to push the mass limits of disk surveys 
towards the least massive objects in star forming regions. Here we present a Spitzer sur- 
vey for disks around IPMOs in the cr Orionis cluster, which harbours the largest IPMO 
population identified thus far. With a well-explored stellar and br own dwarf population 



(jBeiar et al. 



2002c 



2001 



Caballero et al.l 120071 ) and a likely age of 2-4 Myr (jZapatero Osorio et al. 



Sherry et al. II2004I ). a Ori is ideal for probing the disk frequency as a function of mass 



as well as the longevity of IPMO disks. 
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2. Observations and photometry 



As part of the Spitzer program #30395 (PI: A. Scholz), we obtained deep IRAC images 
in the four channels centered at 3.6, 4.5, 5.8, and 8.0 fim for 18 objects in the cr Ori cluster 
with estimated mass es below 20 M.T,ip. AH targets have been identified initia lly in deep 
photometric surveys (IZapatero Osorio et al.l l2000l . l2002bl : ICaballero et al.l 120071 ) With one 
exception, all objects have been confirmed as likely ultra-low mass young cluster members 
with low-resolution spectroscopy; however, given the low signal-to-noise of the spectra, this 
should not be taken as a definite membership confirm ation. In particular, th e membership 
is under debate for SOriTO (see Sect. S]) and SOri47 jMcGovern et al.ll2004l ). See Table □ 
for more information and references for the target properties. 

The objects, covered with 13 IRAC fields, were observed in a 12 position dither pattern 
with the medium scale and 100 sec integration time per position. In total, this gives 1200 sec 
on-source time per target in each channel, a factor 10 deeper (in signal-to- noise ratio) than 



the G TO/IRAC survey of this cluster (program #37, PI: G. Fazio, see iHernandez et al. 



( 120071 )). For the analysis, we made use of the post-BCD mosaics provided by the IRAC 
pipeline version S15.3.0. While parts of some images are affected by saturation effects and 
straylight from bright stars, in the regions where our targets are located the image quality 
is in most cases excellent. 

All objects are detected in channels 1 and 2; two are not detected in channel 3 while 
four are invisible in channel 4. For the object SOri 54, we used the publicly available GTO 
survey data to measure fluxes in channels 1 and 2, because in our deep images the source is 
'drowned' in the emission from a nearby extended object and two nearby point-sources (all 
within 8"). The object SOri 66 is located on a dark row in our IRACl and 2 images; again 
we obtained fluxes for this object from the GTO data. 

Fluxes were measured by aperture photometry with relatively small apertures (3-5 pixels) 
using routines in daophot. For all objects, the aperture is free from visible neighbours. Aper- 
ture corrections were applied, as given in the IRAC data handbook (V3.0). Fluxes were 
converted to magnitudes using the recommmended zeropoints quoted in the data handbook. 
All IRAC magnitudes for the IPMOs are reported in Table [H 

The dominant error source is background uncertainty. To get a realistic estimate, we 
measured the flux in empty regions in our images using the same method as for the targets 
(>50 measurements in IRACl and 2, >80 in IRAC3 and 4). As expected, these values scatter 
around zero; their rms defines the uncertainty in flux measurements due to background noise 
and small-scale fluctuations, and thus gives a conservative estimate on the photometric error 
for faint sources. We obtain 1.8, 2.2, 12.7, 20.3 ;uJy in channels 1-4, respectively. Note that 
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these values do not correspond to the detection limits, which are better defined using the 
peak count rate of the PSF in comparison with the background scatter in the sky annulus. 
All objects for which fiuxes are given in the tables are detected with a signal-to-noise-ratio 
of at least 4. 

The total errors quoted in Tables [1] are comprised of photometric errors and system- 
atics. The latter ones are the combination of uncertainties in calibration (5%) and aper- 
ture correction (1-2%, from the IRAC data handbook). No colour correction was carried 
out; instead we added a small contribution (0-2%) to the error budget. The upper limits 
quoted in Tables [1] correspond to three times the photometric errors. In channels 1 and 
2, our fluxes are in excell e nt ag reement (within 0.2 mag) with the values reported recently 
by IZapatero Osorio et al.l (120071 ) for a subsa mple of our targets. In chan nels 3 and 4, the 
magnitude differences between our work and IZapatero Osorio et al.l (120071 ) typically remain 
within 0.6 mag, still within the la errorbars. We note, however, that in these two bands our 



fluxes are in most cases substantially lower than the values given by IZapatero Osorio et al. 
J2007h . 



3. Disk Frequency in the Planetary Mass Regime 



We search for mid-infrared emission from dusty disks in our sample by comparing the 
measu red IRAC colours with published values for (disk-less) fleld objects in the same Tes 
range (jPatten et al.ll2006l ). The main discriminators to distinguish between the photospheric 
contribution and the disk excess are the fluxes in the IRAC channels 3 and 4; at shorter 
wavelength the contrast between the photosphere and a possible disk is difficult to detect. 
In Fig. [U we plot IRAC colours (11-14 and 11-13) for all targets (except SOri70, see Sect. 
H]). As abscissa in these plots we use the I — J colour, which can be assumed to be purely 
photometric and thus represents a proxy for T^g, which in turn correlates with mass for 
coeval objects. The plot covers a mass range from ~ 20 to ~ 8 Mjup. IPMOs are shown as 
crosses; the solid lines are linear fits to the colours of the field objects and thus delineate the 
photospheric level in these figures. 

As can be seen in this plot, the majority of the IPMOs have IRAC colours indistinguish- 
able from the photospheric values and can thus be assumed to have either no disk or only 
small amounts of dust in the inner disk. Sources that lack excess in our data include SOri54, 



SOri55, SOri56 for which IZapatero Osorio et al.l (120071 ) reported disk detections based on 
the shallow GTO IRAC images. If we use the IRAC3 colour as disk criterion (upper panel), 
we find that four objects - S0ri71, SOriGO, SOri66, SOri58 - show excess at the 3a level, 
i.e. 4 out of 17 or 24±]^f %. Since disk excesses may appear only beyond 5.8 //m, the result 
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from IRAC3 should be considered to be a lower limit. 



Based on the IRAC4 colour (lower panel), we find four sources - SOri 71, SOri J053949.5- 
023130, SOri 60, and SOri 65 - with > 3a colour excesses, which we consider to be primary 
disk detections^ In all four cases the excess increases significantly in IRAC channel 4 com- 
pared with channel 3, indicating rising fiux levels towards longer wavelengths, a clear signa- 
ture of disk emission. With the exception of S Ori 65, these objects have been published pre 



vious ly as disk-bearing very low mass sources (ICaballero et al.l 120071 : IZapatero Osorio et al. 
20071 ). This gives a disk fraction of 4 out of 14 or 29±J3%. Here we do not count the three 



objects with upper limits well-above the photospheric level, for which we cannot decide if 
they have disk excess or not. Disk frequencies derived from IRAC3 and IRAC4 are thus 
consistent; since the value determined from IRAC4 is likely to be more robust, we put more 
emphasis on this result. We note that the disk fraction in our sample might still be some- 
what higher than given here, due to the combin ed effects of photome tric uncertainties and 
contaminating field objects (see the dicussion in lCaballero et al.ll2007l ). 



We now compa re our disk detect i on ra te with previous results for more massive objects 
in a Ori, as given by lHernandez et al.l (120071 ) based on IRAC data: 15% for HAeBe stars, 27% 
for intermediate-mass T Tauri stars, 36% for T Tauri stars, 33% for bro wn dwarfs. The value 



for bro wn dwarfs is in agreement with the disk fraction of 33% derived by lJayawardhana et al. 



(120031 ) from groun d-based L'-band imaging. A higher brown dwarf disk fraction of 47% has 
been published bv lCaballero et al.l J2007I ). For the IMPO range (objects with M < 20 Mj^p), 
we now derive a disk fraction of 29%, which is compatible with the values for T Tauri stars and 
brown dwarfs within the la uncer tainties. We do not see a tre nd to higher disk frequencies 
in the IPMO range, as claimed by lZapatero Osorio et al.l (120071 ). instead the evidence points 
to comparable disk fractions for planetary mass objects, brown dwarfs, and T Tauri stars, 
i.e. over more than two orders of magnitude in object mass (0.008 ... 2 Mq). 

We note that two of the objects with disk excess, SOri 66 and SOri 71, stand out from 
the rest of the sample, as they show excessively strong Ha emission with e quivalent widths 
of ~ 100 and ~ 700A, respectively ( Barrado y Navascues et al.ll200ll . I2OO2I ). This indicates 
that the presence of a dusty disks is likely accompanied by ongoing gas accretion, causing 
intense Ha emission, as observed in T Tauri stars and brown dwarfs. 



^The object SOri 66, albeit having a high 11-14 colour of 1.88 mag, is not counted here, because it has a 
large errorbar in IRAC4. A particular case may be the object SOri 58, which does have excess in IRAC3, 
but appears below th e phot ospheric level in IRAC4. This object has been identified to have a disk by 



Zapatero Osorio et al.l (|2007f ). based on a IRAC4 flux which is significantly higher than our measurement. 



Possible reasons for these ambiguous findings include significant background variations in that area, which 
might hamper accurate photometry. 
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The case of SOri 70 



The faintest target in our sample, SO ri 70, has been reported origir iaUy as a ~ 3 Mjup 
IPMO with tentative spectral type T5.5 (IZapatero Osorio et al.l l2002bl ) . If confirmed, it 
would be the lowest mass free-floating object found to date. Its mas s would put it close 
to or even beyond the predicted opacity limit for fragmentation (see iBonnell et al.l 120071 . 
and references therein ), and thus poses a challenge for star formation theory. However, 
Burgasser et al.l (120041 ) have questioned the cluster membership (and thus the extremely 
low mass) of SOri 70, and argue that its near-infrared spectr um is perfectly consistent with 
a foreground dwarf with spectral type ~T6-7. In response, iMartiru (120041 ) reinterated the 
claim for cluster membership and youth mainly based on the {H — K) colour of SOri 70, 
which is unusually high for a T dwarf, and is interpreted as evidence for low surface gravity. 
To date, the nature of this object remains unclear. Our deep IRAC images now allow a 
re-investigation of the issue. 

SOri 70 is clearly detected in IRAC channels 1 and 2. In addition, it is detected in 
IRAC3, with a ~ 4cr signiflcance (peak countrate over background noise). The IRAC3 flux 
uncertainty has been derived based on the very fact that the object is just detected - if it 
were > 25% fainter than the estimated 17.2 mag, we would not be able to see it. In|2]we 
plot its IR AC colours vs. spect ral type in comparison with measurements for fleld T dwarfs, 
taken from jPatten et al.l (120061 ). SOri 70 is tentatively plotted spectral type T6 (average of 
the two available literature estimates). 

As can be seen in this flgure, the fleld T dwarfs form a clear sequence, whereas SOri 70 
stands out: In both IRAC colours, it appears to have some excess. Compared with the 
mean trend for the fleld T dwarfs, the signiflcance of this excess is ~ 2a in both IRAC2 and 
IRAC3. This excess cannot be accounted for by uncertainties in the spectral type, since the 
3.6 — 5.8 fim colour saturates in the late T dwarf regime and reaches maximum values of 
~ 1.2 - compared with 1.6 for SOri 70. It should also be pointed out that SOri 70 is brighter 
than the more massive L5 dwarf SOri 67 in IRAC2 and 3, again indicating an unusual SED 
for this object. 

There are two poss ible origins for a mid-infrared colour excess in SOri 70: a) As reported 
by lLeggett et al.l (120071 ). gravity affects the near /mid-infrared colours of mid/late T dwarfs 
in the sense that it can produce a signiflcant excess at 2-6/im for low-gravity objects, b) 
Similarly to the L type IPMOs in a Ori, SOri 70 might harbour a dusty disk, which produces 
the IRAC excess. Both possibilities provide additional evidence for youth, and thus the mid- 
infrared excess bolsters the claim that it is the lowest mass free-floating object identifled 
thus far. In case the mid-infrared excess is due to a dusty disk, this can be interpreted as 
evidence for a star-like infancy. While the uncertainties in the IRAC fluxes are too large 
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at this stage for a definitive answer, these new findings certainly motivate further follow-up 
work on this target. 



5. Summary 



We have conducted a disk survey for isolated planetary mass objects, based on deep 
observations with the Spitzer Space Telescope. Our targets comprise the full sample of 
spectroscopically confirmed IPMOs in the a Orionis cluster. As a pioneering study, the 
results presented here are not intended to be the last word on the subject. Instead, we 
intend this Letter to serve as an important step towards a more profound understanding of 
the nature and origin of the lowest mass objects found in isolation. 

A substantial fraction of our targets exhibits mid-infrared excess emission indicative for 
the existence of a dusty disk. From our IRAC data, we derive a disk fraction of 29±J3% 
(4/14) at 8.0 fim. A similar value is found at 5.8 fim. These results are fully consistent 
with the disk frequencies derived for T Tauri stars and brown dwarfs in the same cluster. 
The detection of dusty disks around IPMOs in a Ori is another firm evidence for their 
nature as lowest mass members of this young cluster. The finding clearly establishes that 
objects at the bottom of the IMF can harbour dusty disks with lifetimes of at least 2-4 Myr, 
the most likely age of the a Ori. Thus, our results fit into previous claims for a T Tauri- 



20061 : iJayawardhana fc Ivanovll2006l ). Disk fractions and thus lifetimes are similar for objects 
spanning more two orders of magnitude in mass (0.008 to 2Mq), possibly indicating that 
stars, brown dwarfs, and IPMOs share a common origin. Star formation theory thus has to 
account for a number of objects with masses below the Deuterium burning limit. 

Our sample includes SOri70, the faintest candidate member so far in cr Ori, which is 
detected in our IRAC images from 3.6 to 5.8 fim. Compared with field T dwarfs, the source 
stands out based on the IRAC colours, with excesses at 4.5 and and 5.8 fim at a 2a level. 
This may be an indicatio n for youth and possibly disk occurence for an object with an 
estimated mass of 3Mjup fIZapatero Osorio et al.ll2002bl ). Uncertainties are large, though, 
meaning that this should be seen as motivation for further follow-up rather than a definite 
confirmation of the 'cluster planet' nature of this object. 



like phase in the planetary mass reg ime (INatta et al.ll2002l : iLuhman et al.ll2005l : lAUers et al. 
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Table 1. Targets and Spitzer photometry in the aOri cluster 



Name 


SpT'' 


lb 


I 




3.6 fj. 


4.5 At 


5.8A1 


8.0(1 


comment 






(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




S0ri71 


LO 


20, 


.02 


2, 


.88 


15.24 ±0.06 


14.83 ±0.06 


14.29 ±0.08 


13.84 ±0.13 


disk 


SOri47 


LI. 5 


20, 


.53 


3, 


.15 


14.87 ± 0.06 


14.91 ±0.06 


14.77 ±0.11 


14.90 ±0.33 




SOriSO 


M9 


20, 


.66 


3, 


.12 


15.66 ±0.06 


15.64 ± 0.06 


15.95 ±0.30 


15.63 ±0.70 




SOri51 


M9 


20, 


.71 


3, 


,50 


15.42 ±0.06 


15.30 ±0.06 


15.19 ±0.15 


15.35 ±0.51 




SOri53 


M9 


21, 


.17 


3, 


,28 


16.02 ±0.06 


15.84 ±0.06 


16.69 ±0.63 


15.47 ±0.58 




SOri54 


M9.5 


21, 


.29 


3, 


.30 


15.8 ±0.1'= 


~ 15.6'=''' 


- 15. 3<^ 


> 15.1 




SOri55 


M9 


21, 


.32 


3, 


,10 


16.41 ±0.06 


16.19 ±0.07 


16.18 ±0.37 


15.83 ±0.90 




SOri56 


L0.5 


21, 


.74 


3, 


,30 


16.12 ±0.06 


15.79 ±0.06 


15.62 ±0.22 


15.27 ±0.48 




SOri58 


LO 


21, 


.90 


3, 


.30 


16.43 ±0.06 


16.19 ±0.07 


15.52 ±0.20 


- 16.7'* 


disk? 


SOriJ053949. 5-023130 


NA 


22, 


.04 


3, 


,15 


16.49 ±0.06 


16.32 ±0.07 


15.73 ±0.25 


15.05 ±0.38 


disk 


SOri60 


L2 


22, 


.75 


3, 


,58 


16.61 ±0.06 


16.34 ±0.07 


15.34 ±0.18 


14.34 ± 0.20 


disk 


SOri62 


L2 


23, 


.03 


3, 


.59 


16.50 ±0.06 


16.36 ±0.07 


16.02 ±0.32 


^ 17.2'* 




SOri65 


L3.5 


23, 


.23 


3, 


.33 


16.72 ±0.07 


16.60 ±0.08 


16.05 ±0.33 


15.05 ±0.38 


disk 


SOri66 


L3.5 


23, 


.23 


3, 


,40 


17.41 ± O.l'^ 


16.87 ±0.1'= 


16.05 ±0.33 


15.53 ±0.62 


disk? 


SOri67 


L5 


23, 


.40 


3, 


.49 


17.75 ±0.10 


17.64 ±0.16 


~ 17.6'i 


> 15.1 




SOri68 


L5 


23, 
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Fig. 1. — IRAC col ours for IMPOs in aOri (crosses) in comparison with field M/L dwarfs 
(jPatten et al.ll2006l . solid line). Objects with > Sex excess with respect to the field dwarfs 
are marked with squares. 
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Fig. 2. — The IRAC colours of SOriTO (square) compared with field T dwarfs (jPatten et al, 
20061 . small crosses). 



